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Abstract. We present a survey of the polarization properties of mm-wavelength methanol masers, comprising both 
classes, and transitions from 84.5 to 157.0 GHz. Linear polarization is found in more than half of the observed 
objects, and circular polarization is tentatively detected in two sources. Class I and Class II CH3OH masers show 
similar polarization properties. The largest linear polarization is found in the 133 GHz Class I maser towards L379 
(39.5 %), and in the 157 GHz Class II maser towards G 9.62+0.19 (36.7%). The spectral profiles of the polarization 
angle of Class I masers are mostly flat, except for two sources showing a linear slope. Since the mm-line methanol 
masers are expected to be weakly (or not) saturated, we suggest that the stronger fractional polarizations found 
by us are enhanced by anisotropic pumping and radiative losses. In NGC 7538, we find, for both maser classes, a 
good agreement between our polarization angles, and those measured for the submillimeter dust continuum. This 
can be taken as evidence for magnetic alignment of dust grains. It is also possible that an unsaturated maser with 
equally populated magnetic substates simply amplifies polarized continuum seed radiation. For Class II masers, 
the polarization properties of the various velocity components towards a given source with detectable polarization 
are quite homogeneous. A possible explanation is discussed. Since methanol is non-paramagnetic, the circular 
polarization of the unsaturated maser emission can only be due to variations of the angle between the magnetic 
field and the line of sight along the maser propagation path. 
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1. Introduction 

The polarization of astronomical masers is a complex 
field of current research. Whereas the polarization proper- 
ties of non-amplified radiation simply derive from asym- 
metries of the emitting volume, masers introduce addi- 
tional aspects due to their excitation and propagation. In 
terms of excitation, methanol masers can be classified into 
two well distinguished classes : Class I sources (prototype 
Orion BN/KL) are loosely associated with, but clearly 
offset from ultracompact Hii-regions. They may form in 
outflows and cloud-cloud encounters. These masers are 
thought to be pumped by collisions into high rotational 
levels from where they decay spontaneously into the up- 
per level of the masing transition (Cragg et al., 1992). 
Class II sources (prototype W3 OH) are often coincident 
with known OH masers and Hii-regions, and may thus 
occur in the compressed envelopes exterior to the ionized 
gas. Their masing transitions suggest pumping by FIR 
photons to torsionally excited levels (Sobolev & Deguchi, 



1994). 

Only two studies of methanol maser polarization proper- 
ties are published (Caswell et al., 1993, Koo et al., 1988). 
Both surveys were made at 12 GHz and cover only a small 
sample of Class II masers. Linear polarization of up to 
10% towards 4 (out of 11) and 2 (out of 5) sources re- 
spectively was detected in these two surveys. 

Here we present a survey of polarized millimeter-wave 
methanol masers, including both maser classes. The polar- 
ization of Class I masers was observed for the first time. 
The wealth of strong methanol masers detected at sev- 
eral mm transitions combined with the efficiency of a new 
polarimeter on the 30m telescope were instrumental in im- 
proving upon the two cm surveys with respect to the num- 
ber of detections of polarization, the better polarization 
sensitivity, and the advantage of doing all Stokes param- 
eters simultaneously. 

2. Instrumentation and Observations 
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The observations were done at the IRAM 30m telescope 
on Pico Veleta in February and May 2002. We used the ob- 
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servatory's 3 mm and 2 mm receivers (brought into phase 
with a common local oscillator reference) together with 
a versatile polarimeter working in the intermediate fre- 
quency band near 150 MHz. The inherent frequency agility 
was essential in measuring the polarization of several tran- 
sitions quasi simultaneously. The polarimeter (described 
in more detail in Thum et al., 2003) derives all four Stokes 
parameters from the output signals of the two orthog- 
onally polarized receivers, after insertion of phase lags 
of 0°, 90°, 180° and 270°. For each of these phase lags, 
the coherent superposition of the two signals yields an 
output signal whose power contains two terms, (i) the 
power from the two receivers, and (ii) a correlation term 
(like in an adding interferometer). The correlation terms 
are extracted by considering power differences between 
the four output channels, and yield Stokes I, U and V. 
Stokes Q is derived from the difference of the total pow- 
ers of the two receivers. Small phase offsets between both 
receivers (of order ^ 5°/hour) were regularly (typically 
once per hour) measured with a signal generator, and re- 
moved with a phase shifter acting on the local oscillator 
signal. The remaining instrumental polarization was reg- 
ularly calibrated by means of unpolarized sources (mainly 
planets), for which all the output channels must yield the 
same signal strength. Any residual instrumental polariza- 
tion on the optical axis was thus removed. The remaining 
polarization sidelobes mainly arise from the optics in the 
receiver cabin, and are well known (Thum et al., 2003). 
For the observations presented here, the beam polariza- 
tion only limits the accuracy of the results for mispointed 
observations (see below) , or masers spots scattered across 
a significant fraction of the polarization sidelobes. 

Since the polarimeter's six output signals were analysed 

by the facility autocorrelator, we were restricted to use a 
39 kHz channel spacing, together with a 17.5 MHz band- 
width. We used a parabolic apodization (Welch time lag 
window), yielding a half-power width of the spectral re- 
sponse function of 1.59 x the channel separation, which is a 
good compromise between spectral resolution and sidelobe 
suppression. The resulting velocity resolutions are given in 
Table 1, together with the characteristics of the observed 
maser lines. The observed linewidth was approximately 
deconvolved from the spectral response to the parabolic 
apodization, using a Gaussian fit to the main lobe of the 
latter. 

Our source selection is based on the strongest mm- 
wavclcngth masers found in the available surveys. Class I 
masers at mm wavelengths have been observed at 84.5 
& 95 GHz (Val'tts et al., 1995a) and 133 GHz (Slysh et 
al., 1997), mm Class II masers at 107 GHz (Val'tts et al., 
1995b) and 157GHz (Slysh et al., 1995). Many of the po- 
sitions of Class II masers are from the 6.6 GHz survey 
of Mcnten (1991). Some mm methanol masers have been 
observed at the same velocities, but lower frequencies, 
with better positional accuracy (Kogan & Slysh, 1998, for 
44 GHz Class I masers with the VLA, Minier et al., 2000, 



for 6.7 & 12.2 GHz Class II masers with the EVN and 
VLB A). 

In order to avoid false polarization detections, we retain 

only those masers that are not affected by the instrumen- 
tal polarization sidelobes, applying the following criteria: 

(1) Any mispointing results, even for an unpolarized sig- 
nal, in an instrumental polarization. Since the latter is 
constant across the intermediate frequency band, a scaled 
copy of the Stokes I signal would appear in Stokes Q, U 
and V. In that case, the scaling factors must be compati- 
ble with the instrumental signature, known from our po- 
larization sidelobe measurements. If this is the case, those 
data are discarded. For two or more masers offset within 
the beam, the result is not unique anymore, due to the 
convolution with the polarization sidelobes. 

(2) Wherever interferometric observations of mm- and 
other methanol masers indicate that the maser spots are 
scattered across the polarization sidelobes of the optical 
system of the 30 m telescope, we discard the data. This 
is the case e.g. for the two 84 GHz maser components in 
DR21(0H) (Batrla & Menten, 1988). 

The lincvir polarization given here is 

^ _ V Qline + ^line f-,\ 
Pl^ - J (1) 

-'line 

i.e. the linear polarization is quoted against the spectral 
line Stokes I flux density (after subtraction of a contri- 
bution from the continuum emission). The relative error 
of the fractional linear polarization is (to first order in 
Pl) given by the signal-to-noise ratio of the Stokes I mea- 
surement. Continuum polarization measurements, as com- 
pared to spectral line polarimetry, are thus more precise, 
owing to the larger bandwidth. A few sources display con- 
tinuum emission sufficiently strong to contribute to the 
signal in our quite narrow spectral band. The importance 
of dust polarization at mm wavelengths is still a matter 
of debate. Linear polarization (up to 8 %) of the dust 
continuum at A 1.3 mm was observed towards the Orion 
cloud core, offset from the BN/KL position (Leach et al., 
1991). Towards BN/KL, the fractional polarization drops 
to 2.7% ("polarization hole"). In our case, a mispoint- 
ing of 10" can also lead to a linear polarization of up to 
8 % (depending on the direction of the mispointing) . The 
fractional continuum polarization may also be affected by 
the measurement inaccuracy of Stokes I due to fluctua- 
tions of the atmospheric emission (the IF polarimeter is 
not designed to measure weak continuum polarization lev- 
els). Stokes Q, U and V are not concerned, since they 
result from differential and correlation measurements, re- 
spectively, so our continuum polarization angle measure- 
ment remains reliable - on condition that the continuum 
emission is unresolved. 

The statistical bias, introduced in the determination 
of ^Q-^ + was subtracted from the percentage polar- 
ization (see Wardle & Kronberg, 1974). The standard de- 
viations of the linear polarization and its position angle 
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Table 1. The observed methanol masers. All masers are in the torsional ground state. Rest frequencies are from Anderson et 
al. (1990). 



Ju 


Ku 




Jl 


Ki Symm. 


Frctjucncy 


En 


El 


Vcl. resolution 


Maser 










Substate 


[MHz] 


[cm-1] 


[cm-i] 


[kms^i] 


class 


5 


-1 




4 


E 


84521.210 


27.417 


24.597 


0.221 


(I) 


8 







7 


1 A+ 


95169.440 


14.905 


11.705 


0.196 


(I) 


3 


1 




4 


A+ 


107013.850 


19.704 
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0.140 


(I) 


6 







6 


-1 E 


157048.620 


42.331 


37.092 


0.119 


{II) 



Table 2. Detections (> 3arms) of linear polarization towards Class I CH3OH masers. Only velocity channels with the most 
significant polarization arc given here (with velocities in column 10). For the other velocities, see Figs. 2 to 6). The maser 
polarization is given in columns 8 (fractional linear polarization pl) and 9 (polarization angle Xj East from North). The maser 
fluxes and half-maximum full widths are derived from least-square fits of Gaussian components; the resulting linewidths are 
corrected for the instrumental broadening (for details see text). 



Source 


Ref. 


a(2ooo) 


•^(2000) 


Frequency 




Avfwhm 


PL 


X 


VlBT 










[GHz] 


[Jy] 


[kms-i] 


[%] 


n 


[kms-^] 


OMC-2 


(1) 


05''35™2r!07 


-05°09' 52'.'5 


132.891 


37.7 


0.36 


8.94 ± 0.91 


22.1 + 2.6 


11.21 


S231 


(2) 


05 39 13.06 


+35 45 51.3 


95.169 


12.7 


0.51 


3.79 ± 0.58 


118.1 + 4.5 


-16.82 


NGC 2264 


(1) 


06 41 08.02 


-h09 29 40.0 


132.891 


14.0 


0.62 


10.6 ±1.7 


106.4 + 4.1 


7.22 


M8E 


(2) 


18 04 53.76 


-24 26 35.3 


132.891 


24.6 


0.63 


33.2 ± 1.2 


154.56 + 0.8 


10.93 


W 33- Met 


(1) 


18 14 11.08 


-17 55 57.4 


132.891 


8.8 


0.79 


22.8 ± 3.0 


154.5 + 3.4 


32.95 


L379 


(1) 


18 29 24.69 


-15 15 19.0 


132.891 


3.8 


0.36 


39.5 + 6.0 


133.7 + 3.4 


19.18 


W51-Met2 


(3) 


19 23 46.50 


+14 29 41.0 


84.521 


6.0 


0.38 


9.2 ± 2.0 


23.6 + 5.6 


56.12 










132.891 


3.0 


0.49 


20.1 ± 1.1 


10.2 + 3.4 


53.28 


DR21-W 


(1) 


20 38 54.72 


-h42 19 22.4 


84.521 


85.6 


0.35 


2.77 + 0.15 


48.1 ± 1.4 


-3.15 










132.891 


72.5 


0.42 


4.26 + 0.50 


7.7 + 3.0 


-2.44 


W75 S(3) 


(3) 


20 39 03.50 


+42 25 53.0 


84.521 


14.3 


0.28 


9.4 + 2.0 


7.8 + 5.7 


-3.62 


NGC 7538- Metl 


(3) 


23 13 46.40 


+61 27 33.0 


84.521 


3.8 


0.70 


5.5 + 1.6 


161.2 + 5.0 


-58.00 










95.169 


14.1 


0.82 


14.5 ± 1.3 


7.2 + 2.6 


-57.42 










132.891 


3.5 


0.81 


24.6 ± 1.8 


53.3 + 1.9 


-57.35 


References: (1) VLA positions of the 44 GHz maser, Koj 


;an & Slysh (1998), (2) VLBI positions of the 6.7 & 12.2 GHz 


masers. 


Minier et al. (2000), (3) positions of the 95 GHz survey of Val'tts (1995a). 










Table 3. As Table 2, but for Class II CH3OH masers. The distinct main velocity components are listed separately. 




Source 


Ref. 
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'5(2000) 
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X 
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[GHz] 


[Jy] 


[kms-i] 


[%] 


n 


[kms-1] 


NGC 7538-IRSl 


(1) 


23''13™45!36 


+61°28'09'.'7 


107.014 


< 16.3 


< 2.34 


3.63 ± 0.97 


104.5 + 6.8 


-56.30 










107.014 


< 11.8 


< 1.57 


5.7 + 1.2 


116.2 + 5.8 


-58.59 










107.014 


8.9 


0.88 


6.1 + 1.6 


112.6 + 6.7 


-60.56 


CepA 


(1) 


22 56 18.10 


+62 01 49.5 


157.049 


19.5 


0.48 


12.3 + 1.7 


86.4 + 4.0 


-3.84 










157.049 


40.6 


0.34 


14.10 + 0.94 


87.2 ± 2.0 


-2.80 










157.049 


17.0 


0.57 


14.0 + 2.0 


79.0 ± 4.2 


-1.75 










157.049 


10.2 


0.58 


10.6 + 2.7 


81.2 + 6.9 


-1.30 


G9.62-h0.19 


(1) 


18 06 14.80 


-20 31 32.0 


157.049 


6.4 


0.88 


31.5 + 3.7 


142.9 ± 3.1 


1.55 










157.049 


3.5 


0.66 


36.7 + 5.7 


135.1 + 4.8 


-0.09 



Note: The masers of NGC7538-IRS1 are blended with thermal emission, only upper limits can be given for flux and linewidth. 
References: (1) VLBI positions of Minier et al. (2000) 
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were derived from the polarimeter's output signals by er- 
ror propagation. 

3. Results 

In the foUowing, the detections of methanol maser polar- 
ization are discussed within the context of the environ- 
ment of their star-forming regions. The two maser classes 
are well separated, i.e. each of the transitions summa- 
rized in Table 1 show maser action only in one of the 
two CH3OH maser classes. We find (all transitions con- 
founded) polarization in 10 out of 14 Class I sources sat- 
isfying our selection criteria, and in 3 out of 7 Class II 
sources. We cannot break down this statistics into spe- 
cific transitions, since not all sources were observed at 
all frequencies of Table 1. The distribution of polariza- 
tion detections is shown in Figure 1. The detections of 
linear polarization range between 2.8% and 39.5%. The 
velocity profiles of polarization angles are either flat, or 
show a linear slope. Circular polarization is tentatively 
3(7rms) detected towards two Class I masers. Our re- 
sults are summarized in Table 2 for Class I mm CH3OH 
masers, respectively in Table 3 for their Class II coun- 
terparts. Given the non- uniformity of our data set, and 
the fact that the accuracy of the polarization measure- 
ments depends on the signal-to-noise ratio, the upper lim- 
its for the non-detections depend much on the frequency 
and the weather conditions. Typical upper limits for non- 
detections range from = 0.5 % to about 10 %. 

3.1. Class I Masers 

The polarization spectra for Class I masers are shown in 
Figures 2 to 4, for 84 GHz, 95 GHz, and 133 GHz, respec- 
tively. The spectra of the fractional linear polarization (in 
percent) and its position angle are for the velocity range 
indicated by the grey-shadowed area in the Stokes I spec- 
tra. All polarization angles are given in counterclockwise 
sense (E from N). 

3.1.1. OMC-2 

We observed this maser, located in a cluster of in- 
frared sources and far-infrared condensations with mul- 
tiple outflow activity, in the 133 GHz 61 5o tran- 
sition of E-type methanol (flrst detection by Slysh et 
al., 1997). Other Class I transitions were previously ob- 
served, but with unreliable astrometry. The maser peaks 
at 11.35kms~^ , the strongest polarization (pl = 8.9%) 
is observed at 11.21 kms^^ . At flrst sight, the profile of 
the linear fractional polarization looks like a scaled copy 
of the Stokes I line profile. However, the polarized fea- 
ture is significantly narrower than the maser component 
(0.37kms~^ vs. 0.69 kms~^ , respectively). An instrumen- 
tal origin can also be ruled out, since criterion (1) of 
section 2 does not apply. The polarization angle (22°) 
remains constant across the line profile. As a matter of 
fact, the outflow driven by the source FIRS and mapped 
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Fig. 1. The distribution of detections of lineax polaxization 
towards Class I masers (top, for a total of ten sources) and 
Class II masers (bottom, for three sources), binned into 5% 
intervals. 



by Williams et al. (2003) has a position angle of 30°; 
the brightest maser spots of the 7o 61 ^4+ CH3OH 
maser (at 44 GHz, mapped with the VLA by Kogan & 
Slysh, 1998) are on a straight line of position angle 32° (E 
from N). Both position angles are close to our polariza- 
tion angle. This interesting finding, schematically shown 
in Figure 7, will be discussed below. The position and 
velocity of the maser spike corresponds to one of the com- 
ponents of the 44 GHz maser (component number 7 of 
Kogan & Slysh). The thermal linewings make it impossi- 
ble to know whether we also get weaker maser emission 
at the velocities of the other VLA maser spots, since the 
velocity dispersion of the latter is merely 0.23 km s~^ . If 
so, they could produce an important instrumental polar- 
ization, because they are spread over a N-S line extending 
over 15" , i.e. as large as the telescope's main beam at that 
frequency. However, this would not be compatible with our 
flat polarization angle profile, since the sidelobe polariza- 
tion pattern is significantly different in Stokes Q and U 
(see Thum et al, 2003). 
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Fig. 3. Class I methanol masers at 95 GHz. 



20 



10 - 



— 



50 - 



1 '"''il 

N H 



-50 



-59 



-58 -57 
V|3r [km/s] 



-56 



3.1.2. S231 



The 95 GHz 8o li A+ CH3OH emission, first observed 
by Val'tts (1995a), is at -16.8 kms-^ , with a 3.8% linear 
polarization. Its polarization angle across the main com- 
ponent shows a clearly defined gradient, a finding that will 
be discussed in the next section. A second, weaker maser 
component, not or weakly polarized, is at — 15.2kms~^ , at 
the blue-shifted end of the cluster of 6.7 GHz 5i — > 60 ^4+ 
CH3OH masers (EVN observations of Minier et al., 2000). 



3.1.3. NGC2264 

The 133 GHz 6_i -boE CH3OH maser, first observed by 
Slysh et al. (1997), is located in a core of this massive star- 
forming region, that is associated with the bright infrared 
source IRS 1 (Allen, 1972). The position angle of the lin- 
ear polarization (^ 110°) remains remarkably constant. 
Schreyer et al. (1997) decomposed the outflow activity of 
IRS 1 into two flows. They infer that one of them, extend- 
ing in North-South direction, is driven by IRS 1 and seen 
nearly pole-on. The 133 GHz CH3OH maser is located in 
that outflow. 



3.1.4. M8E 



We find in the 133 GHz line linear polarization as high 
as 33%. Its position angle of 155° is flat across the spec- 
tral profile. M8E harbours a strong 44 GHz maser (Kogan 
& Slysh (1998), whose two components are 1'.'7 apart 
(2500 AU at a distance of 1.5 kpc, see Val'tts, 1999a). The 
133 GHz maser is at the radial velocity of the strongest 
44 GHz component {v\si = 10. 94 km s^^, the weaker one 
at v\sr — 10.45 kms^^ does not appear at 133 GHz). At 
10. 84 km s~^, we flnd a circular polarization of pc = 
(—7.1 ± 1.0) %, which cannot be due to the sidelobe po- 
larization alone (Fig. 10). The maser is most likely not 
associated with the CO outflow (position angle 54°, W 
from N, Mitchell et al., 1992), since the velocities of the 
quiescent gas and the methanol maser are close to each 
other (since the outflow almost appears pole-on, a large 
velocity offset would be expected for an associated maser). 
Consequently, Val'tts (1999a) suggests that the Class I 
maser is located in the front of the bipolar outflow at an 
interface with the surrounding dense molecular gas. 



H. Wiesemeyer et al.: The polarization of mm methanol masers 



7 




H. Wiesemeyer et al.: The polarization of mm methanol masers 




20 "■■■Z5 
V|sr [km/s] 




'50 "■■60 
V|sr [km/s] 



60 - 



40 



20 



50 



-50 



18 



19 20 
V|3r [km/s] 



21 



20 



10 - 



50 



-50 



52 54 56 

V|3r [km/s] 



100 - 




-5 0" 
V|sr [km/s] 




-60 -55 
Visr [km/s] 



50 



-50 



-4 



4' 4' 4' 



4- 4- 4' 



-3 -2 
V|3r [km/s] 



-1 



30 



20 



10 - 



50 



- 



-50 - 



-60 



-58 -56 
V|,, [km/s] 



H. Wiesemeyer et al.: The polarization of mm methanol masers 



9 



3.1.5. W33-Met 

We measure, for this 133 GHz maser, a fractional linear po- 
larization of 23 % at 32.9 km s~^ . The radial velocity coin- 
cides with that of the southern component of the two VLA 
44 GHz maser spots of Kogan & Slysh (1998). Another po- 
larized feature, with ~ 10 %, is at a radial velocity of 
about 37kms~^ , where the 44 GHz maser has a second 
component, at about 32" northwards. This feature coin- 
cides with a broader CH3OH emission, of equal flux den- 
sity. Val'tts (1999b) suggests that this (spatially extended) 
component is a quasi-maser with a low gain and weak line 
narrowing. Thus, the linear polarization measured by us 
(although barely significant, ^ 3(7j.,ns) towards this com- 
ponent almost certainly has an important instrumental 
contribution. 

No spatially resolved map of the outflow exists. Val'tts 
(1999b) find the Class 1 maser to be associated with a CS 
condensation, containing a very young protostar. 

3.1.6. L379 

The 133 GHz maser is located 30" to the north of IRAS 
source L379-IRS1. From the absence of linear polariza- 
tion in the thermal red-shifted linewing, it is excluded that 
the measured linear polarization of 39 % at a radial veloc- 
ity of 19.2 km s~^ is due to the telescope's sidclobc polar- 
ization. The polarization angle has a flat spectral profile 
with X = 134°. This 133 GHz maser concides with the 
northern A800/xm continuum peak, discovered by Kelly 
& MacDonald (1996), and is probably associated with 
their CO outflow, whose direction is at 75°. It is likely 
that the red-shifted linewing consists of a blend of (quasi- 
)maser emission. As a matter of fact, we find an even 
stronger fractional linear polarization with a similar po- 
larization angle at about 19.6 km s~^ , but of lower signifi- 
cance (3-4 CTrms)- The 44 GHz CH3OH emission (Kogan & 
Slysh, 1998) shows a total of 11 masers, scattered across a 
region of 46" in diameter, in two groups at radial velocities 
of 18.8 -20.8 km s"^ and 16.3 - 18.8 kms'^ respectively. 
From the latter group, we detect CH3OH emission, whose 
line profile can be decomposed into two Gaussian spectral 
components, a large one (linewidth 7.0kms~^ FWHM), and 
a narrower one, with a linewidth of 2.1kms^^(at best 
a quasi-maser without substantial line narrowing). The 
polarization of these components, most likely scattered 
across our beam, is strong across the whole line profile. 
Stokes U is an inverted and scaled copy of Stokes 1, and 
the signs of the other Stokes parameters agree with the 
instrumental signature. We thus discard this polarization 
signal as of instrumental origin. 

3.1.7. W51-Met2 

The 84 GHz line of this CH3OH maser is blended with a 
broad thermal emission. However, the linear polarization 
disappears across a large part of the thermal profile, ex- 
cept for the single-peaked maser line at 56.1 km s~^ , and 



another component at 53.3 km s~^ (possible maser emis- 
sion masked by the thermal emission, cf. Fig. 2). The 
spectral profile of the polarization angle of the 84 GHz 
maser spike has a linear slope, as in the case of S231. 
Jumps in the polarization angle of 90° at certain ve- 
locities arc due to a statistical bias in the weak-signal 
limit. The thermal emission of the 133 GHz maser is 
strongly polarized (between 52.0 and 55.2 km s~^ , reach- 
ing a maximum of 20.1% at 53.3 km s~^ , with a polar- 
ization angle of —5° to 30° in this velocity range). The 
line peak can be decomposed into a broad component 
(linewidth 4.5 km s"-'^ fwhm) and a narrow one (linewidth 
0.482 km s~^ ). We consider it quite possible that this ex- 
traordinary, broad polarization profile may be due to an 
anisotropic optical depth (Goldreich & Kylafis, 1982). The 
amplification of polarized background continuum serving 
as "seed" is also possible, given the fiux density ratio of 
-Fiino/^cont = 1 1 at 53.3 km s^^ (or 2.8 if we only consider 
the narrow component). However, the uncertain astrom- 
etry and extent of the thermal emission could introduce 
an important instrumental contribution, too (as discussed 
in the previous section), which makes impossible an un- 
biased detection of the polarization of the 133 GHz maser 
at 56.6 km s^^ (the linear polarization is not significant, as 
opposed to the 84 GHz maser). 

As an indication of the possibility of polarized seed 
radiation, wc detect a fractional linear continuum po- 
larization of PL = (4.7 ± 1.2)% at 84 GHz, and pl = 
(13.7±2.1) % at 133 GHz (with x = 18?4, not too far from 
the submm measurement, Chrysostomou et al., 2002), re- 
spectively. Note that the errors in pl do not contain the 
uncertainty due to atmospheric fluctuations (as discussed 
in Section 2). 

The host source is the W51f subregion of the W51 A 
star forming cloud (sec Sievers et al., 1991, for a A870/im 
and A 1.3 mm map). Submillimeter continuum polarime- 
try (Chrysostomou et al., 2002) reveals a marginal linear 
dust polarization of 2.5%, at a position angle 24°. The 
authors conclude that the overall magnetic field is rather 
weak, and becomes only important close to the sites of star 
formation. The lack of evidence for a continuum signal at 
longer wavelengths made it difficult to reveal the nature 
of this maser's environment. Hodapp & Davis (2002) re- 
cently found a spatially unresolved, but strong H2 (1 — 0) 
S(l) line (A 2.12 /xm) about 10" eastwards of the methanol 
maser position. Taking into account the uncertain astrom- 
etry (no interferometer observations are available), this 
emission could come from the same shocked gas that har- 
bours the methanol maser, in the immediate vicinity of a 
star enshrouded by the molecular and dust cloud. 

3.1.8. DR21-W 

DR21-W is one of three sources in which we heave mea- 
sured and detected linear polarization in both the 84 GHz 
and 133 GHz line. They both have flat polarization an- 
gle profiles (with position angles 48? 1 and 7? 7 respec- 
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tively). They are single-peaked, and are rather weakly 
polarized (about 3 — 4%). At 133 GHz, we find a sig- 
nificant circular polarization of pc = (3.52 ± 0.68) % at 
Vist: = — 2.63knis~^ (Fig. 10), slightly offset from the ve- 
locity of maximal linear polarization (— 2.44kms~^ ). As 
in the case of M8E, this detection can not be entirely at- 
tributed to the sidelobe polarization. 

Kogan & Slysh (1998) observed the 7o - 6iA+ 
methanol mascr at 44 GHz with the VLA, and find two 
major components separated by 2'.' 8, whose radial veloc- 
ities are close to ours. The environment of this maser is 
fairly well known: Davis & Smith (1996) obtained high- 
resolution images in H2(2.12/xm emission). Their spec- 
trum closest to the methanol maser is best described by 
a C-typc planar shock. The submillimctrc continuum po- 
larimetry (Minchin & Murray, 1994) reveals a linear mag- 
netic field, in position angle (114°) intermediate between 
the cloud core and the dominant outflow axis (70°, Garden 
et al., 1991). This polarization angle, however, does not 
agree well with that of either of the two maser lines. Liechti 
& Walmsley (1997) observed the thermal 2k Ik emis- 
sion of A+ and £^-type methanol. The methanol maser 
shown here coincides in position with a thermal clump 
in 2i — > E emission, whose major axis is at position 
angle 45° - remarkably close to our 84 GHz polarization 
angle. The thermal methanol emission also appears in 
our 84 GHz and 133 GHz spectra as a broad, unpolarized 
pedestal underneath the maser lines. 

3.1.9. W75 S(3) 

In this source, located in the same ridge of molecular gas 
as the DR21 region (Wilson & Mauersberger, 1990), the 
84 GHz masers consist of a narrow main component, and 
a blend of several, weaker (quasi-) maser components, to- 
gether with a broad underlying thermal emission. Only 
the polarization of one of the weaker components is sig- 
nificant, the other maxima of the polarization profile are 
far less significant (see Fig. 2). It is difficult to distinguish 
this feature from an intrinsic or instrumental polarization 
of the underlying, broad emission. The profile of the polar- 
ization angle oscillates between 50° and 80°. At 133 GHz, 
no significant linear polarization was detected. 

3.1.10. NGC7538-Metl 

This is the only source in our sample with a bona-fide de- 
tection of maser polarization in each of the three Class I 
maser frequencies. The polarization angles at the maser 
peaks are different, their velocities are within 0.8kms~^ . 
This difference is most likely intrinsic, since the polariza- 
tion angle profiles themselves have a peak-to-peak varia- 
tion of up to 20°. The strongest polarization is found at 
133 GHz (24.6%), resulting in the most significant polar- 
ization angle measurement among the three transitions. 
The polarization angle has a remarkably flat spectral 
profile (position angle 53?3 at the 133 GHz maser spike. 
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Fig. 5. The 107 GHz polarization spectra of NGC 7538-IRSl, 
a Class II CH3OH maser. 

as compared to 161?2 at 84 GHz, and 7?2 at 95 GHz). 
Although part of the polarization could come from the 
underlying broad emission (at 84 GHz and 133 GHz), it 
is unlikely that the latter reaches a fractional polariza- 
tion as high as 25 %. The A 870 /xm continuum polarization 
is from magnetically aligned dust grains (Momosc et al., 
2001). Its position angle (41?8) differs by 11?5 from our 
measurement at 133 GHz (Figure 8). 

3.2. Class II Masers 

The polarization spectra for Class II masers are shown in 
Figures 5 and 6, for 107 GHz and 157 GHz, respectively. 

3.2.1. NGC7538-IRS1 

The 107 GHz maser from this source is located 37" north- 
wards of its Class I counterpart, and coincides with 
the infrared source IRSl. The 107 GHz methanol maser 
has three components, each of similar polarization. The 
107 GHz continuum of the source, as seen in our spectra, 
is not polarized (upper limit pl < 0.4 %). The velocities of 
the maser spikes correspond to those of the 6.7 GHz maser 
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Fig. 6. Class II methanol masers at 157.0 GHz. 

components A,E, and D (or C) of Minier et al. (2000), 
which are 0'.'6 apart (i.e. 1620AU at 2.7kpc distance). 
The polarization angle of the three velocity components 
of the 107 GHz mascr (Figure 5) is 105° — 116°, not too far 
from the position angle of the aligned 6.7 GHz masers of 
subgroup A, and close to the position angle of the equator 
of the cdgc-on circumstcllar disk proposed by Minier et 
al. They suggest that the maser resides at a single radius 
of this disk. This conjecture is corroborated by the detec- 
tion of an elongated structure seen in an intcrfcromctric 
^■^00(1 — 0) map (Scoville et al., 1986), perpendicular 
to the bipolar HE region, and at an angle of 55° to the 
major outflow axis. As the Class I mascr to the south of 
IRS 1, our Class II maser polarization angle agrees more 
or less (within 10° to 20°) with that of the A870/im dust 
emission (Momosc ct al., 2001, see our Figure 8). The lin- 
ear polarization at 107 GHz yields significant results for 
three velocity components (pl = 3.6 to 6.1%, Fig. 5). 
However, given that this fractional polarization is rather 
weak (for the 157 GHZ maser, no significant polarization 
is detected), it cannot be excluded that the whole line 
profile, with both (quasi-) mascr and thermal emission is 
polarized, due to the amplification of a polarized contin- 
uum source. As a matter of fact, the continuum emission 
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at 157 GHz has a linear polarization of pl = (8.4 ± 1.0) %. 
Its position angle (-84?0 ± 2? 7) is close to that of the 
A870/im polarization (— 87?2, Momose et al.), which dif- 
fers from the polarization angle of the 107 GHz mascr by 
about 23°. For the ratio of the line to continuum bright- 
ness temperature we estimate values of 6 for the 107 GHz 
maser, and of 5 for the 157 GHz maser, respectively. 

3.2.2. CepA 

We observed the 107 GHz and 157 GHz masers towards 
CepA East, detecting a significant linear polarization 
only from the latter. This source contains masers at four 
well separated velocities, which agree with the compo- 
nents of the 12.2 GHz maser (VLBI observation of IMinier 
et al., 2000). The fractional polarization of the compo- 
nents of the 157 GHz masers ranges from pl = 10.6% to 
14.1%, and the polarization angle from 79° to 86° (see 
Fig. 11). The significant polarization measurements (i.e. 
PL > 3fT) are clearly correlated with the maser spikes 
(there is stronger, but less significant polarization between 
the spikes). The velocity of the first component listed in 
Tables agrees well with that of component B in the in- 
terferometer map of the 107 GHz maser (Mehringer et 
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al., 1997). Our last two components correspond in both 
velocity and flux density to component A of Mehringer 
et al. (after smoothing to their spectral charactcrsistics), 
which is spectrally unresolved by them, but asymmetric 
due to the blend of the two spikes. Although Mehringer 
et al. infer from their results that the 107 GHz masers 
arise from regions with slightly different physical charac- 
teristics, their polarization properties are rather homo- 
geneous - although the projected distance between com- 
ponents A and B is 0'.'7 (corresponding to 500 AU at a 
distance of 730 pc, Torrelles et al., 1996). The fourth com- 
ponent at visr = — 2.8kms~"'^ does not have a counterpart 
at 107 GHz. At 157 GHz, it is the strongest one (both 
in intensity and polarization). The velocities of all these 
masers are rcdshifted with respect to those of the ambient 
molecular cloud and the expanding bubble of H2O masers 
(Torrelles et al., 2001). According to Minier et al. (2000), 
the methanol masers and the water maser are located in 
the same environment. However, given the complexity of 
the region, it cannot be ruled out that the CH3OH maser 
rather emerges close to the thermal radio jet of Torrelles 
et al. (1996). Codella et al. (2003) found evidence of a well 
coUimated outflow towards the south; the methanol maser 
could be located just at its onset. 

3.2.3. G 9.62+0.19 

This maser has two components on top of a broad ther- 
mal and possibly quasi-maser contribution. As in the case 
of CepA and NGC 7538-IRSl, both velocity components 
of this maser have the same polarization characteristics 
(weighted meanpL = 33.0%, at x = — 140?6). The masers' 
position and velocity coincide with that of a cluster (0'.'5 
in extent) of 6.7 GHz maser emission (Norris et al., 1993). 
The associated UCHii region of this methanol maser is 
continuum source E (Garay et al., 1993), located near 
the edge of a compact, dense and warm molecular cloud 
core (Hofner et al., 1996). Hot thermal methanol emission 
of the 153 14:4 A~ hne was imaged with the Plateau 
de Bure interferometer by Hofner et al. (2001). It coin- 
cides with the continuum source E and the maser position. 
Although the thermal methanol emission is point-like in 
the interferometer map, there may be an extended com- 
ponent missed by the interferometer. As a matter of fact, 
the whole line profile at 157 GHz is polarized, which may 
be due to the sidelobe polarization of a component that is 
resolved in our beam. 

4. Discussion 

Most models of maser polarization to date have been for 
saturated masers (Goldreich et al., 1973, Watson & Wyld, 
2001). The polarization in such models is due to unequal 
populations in Zeeman-split substates of the maser levels 
and one finds that a high degree of maser saturation as well 
as the presence of a magnetic field is a necessary condition 
for linear polarization in this type of model. One can also 
have polarization in other situations, where the radiation 



field pumping the maser or the optical depth in FIR tran- 
sitions connected to one of the maser levels is anisotropic. 
Western & Watson (1983, 1984) investigated the role of 
anisotropic pumping in the polarization of maser lines, 
both for J = 1 — and J = 2 — 1 transitions. Elitzur 
(1993) considered the exponential maser growth in the 
unsaturated phase, where the polarization is produced. 
Subsequently, it has been found (Elitzur, 1996) that for 
non-paramagnetic molecules, where lincwidths exceed the 
Zeeman splitting, the higher-order terms of the Zeeman 
splitting, although small, cannot be neglected, but have 
to be taken into account in order to rule out unphysical 
solutions. A summary of this work can be found in Elitzur 
(2002). 

In the case of the masers observed by us, we feel that 
one should presently consider all options open and concen- 
trate on obtaining better observational constraints on the 
physic;al conditions in the regions where methanol maser 
polarization is observed. Below we comment briefly both 
on what is known about the physical conditions in the 
clumps giving rise to methanol masers as well as on some 
of the current models. 

4.1. Physical conditions in the regions giving rise to 
methanol maser emission 

Studies of Class I methanol masers to date have been ham- 
pered by poorly known collisional rates. However, a re- 
cent study by Pottage, Flower & Davis (2002, 2004) of 
He-CHsOH collisions has improved the situation, and it 
appears as if the estimates used in previous studies (e.g. 
Gragg et al., 1992, Johnston et al., 1992) were broadly cor- 
rect. Calculations by Leurini et al. (2004) using the new 
rates confirm this. 

On the basis of the earlier rate estimates, Johnston et 
al. found, using VLA observations of the 25 GHz tran- 
sitions towards Orion, that the H2 density in the maser 
rc^o'ion should he bcrtween 2 x 10*" and 10^ cm^"* and the 
methanol cohmm density in the range 2 x lO^*' to 2 x 10^^ 
cm~^. The temperature is expected to be of order 100 K 
and the inferred methanol abundance CH3OH/H2 10^^ 
to 10^^. These results are similar to those obtained for 
other CH3OH lines in other Class I sources (e.g. Liechti & 
Walmsley, 1997, in the case of DR21) and one concludes 
that high density (above 10'"' cm^"'), high temperature (of 
order 100 K), and high methanol abundance are charac- 
teristics of Glass I methanol emission. It is notable more- 
over that the estimated abundances are similar to those 
found by Bachiller et al. (1995) towards outflow regions, 
but much higher than those in cold dark clouds (Friberg 
et al., 1988). One interpretation of this is that the ele- 
vated methanol abundance is caused by sputtering of ice 
mantles of dust grains containing methanol in shocks as- 
sociated with the outflows. If so, the alignment due to the 
shock likely has consequences for the magnetic field ori- 
entation and for the escape probability in optically thick 
CH3OH lines. 
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Class II methanol masers (in particular the 6.7 and 
12.2 GHz transitions) have been studied by Sobolev et al. 
(1997a, 1997b) and by Cragg, Sobolev & Godfrey (2002). 
The general conclusion is that the masers form in high 
density regions (10^ to 10^ cm~^) under conditions in 
which the gas temperature is lower than the dust tem- 
perature (both being of order 100 K). There is evidence 
for association with OH maser regions having milligauss 
magnetic fields. There is evidence also that in W30H, the 
maser regions are associated with absorption components 
in non-masing transitions of OH (Guilloteau et al., 1985) 
and methanol (Menten et al., 1986) whose temperatures 
are of the order 100 - 200 K. In both cases (CH3OH and 
OH), one requires high abundances in the range 10"^ to 
10~^ whose explanation (Hartquist et al., 1995) is thought 
to lie in the evaporation of CH3OH from the surface of 
grains heated by the passage of weak shocks. As a matter 
of fact, ISO observations (Gibb et al., 2000) have shown 
that CH3OH is an important constituent of interstellar 
ices associated with protostellar objects such as W33A. 
For low mass stars, Pontoppidan et al. (2003) detect abun- 
dant sohd methanol, too. 

4.2. Models of polarization of saturated masers 

Goldreich et al. (1973) showed that the Zeeman splitting 
Ob in polarized masers has to exceed \/RT, where R is 
the stimulated emission rate of the masing transition, and 
r is the decay rate out of the maser levels. However, this 
result only holds for J = 1 — transitions. For higher 
angular momenta, Dcguchi & Watson (1990) showed that 
this condition has to be replaced by the stronger require- 
ment Ob > R- It is difficult to get even a rough esti- 
mate of R from observable quantities, since the spot size 
of mm CH3OH masers is not well known. Lonsdale et al. 
(1998) deduce from mm VLBI measurements a lower limit 
of several milliarcseconds. In the following, we assume that 
i? ^ r, as for unsaturated maser emission. We estimate for 
the gas density in CH3OH masers values of 10^— lO'' cm~^ , 
where magnetic fields of around a 1 mG are expected. The 
Zeeman splitting is thus about 0.3 Hz. For the Lande fac- 
tor of CH3OH, a value of 0.078 has been used (Jen, 1951), 
or equivalently 4.2 x 10~^ if the Bohr magneton is used 
instead of the nuclear magneton. The loss rate F, on the 
other hand, is around 10~^ — 10~^ s~^, for a collision rate 
of cm~'^nH2, and the above gas density. The Zeeman 
splitting thus exceeds R by several orders of magnitude, 
except if the maser is strongly saturated (which is not the 
case here). 

According to Goldreich et al. (1973), the polarization 
angle is either perpendicular to the projection of the mag- 
netic field on the sky (hereafter B^ky), if sin^ O > 2/3, 
where O is the angle between the magnetic field direction 
and the maser propagation direction, or parallel otherwise. 
However, this result only holds in the limiting case of ex- 
treme saturation T « R and magnetic fields satisfying 
R « Ob. 



Extensions of this model were studied by Nedoluha & 
Watson (1990, hereafter NW90), including lower satura- 
tion and larger angular momenta. Since they are relevant 
for our results, we discuss them in the next section. 

4.3. Extensions of the idealizing maser polarization 
models 

NW90 (their Fig. 3a, b) showed that for all saturation pa- 
rameters R/r and for a maser propagation with sin^ <d > 
2/3 (yielding the largest linear polarization) the polariza- 
tion angle remains (except for a small deflection) roughly 
perpendicular to Bgky. For smaller G, all orientations of 
the polarization angle, from parallel to Bgky for weak satu- 
ration to perpendicular to Bsky for strong saturation, are 
possible. For a J = 2 — 1 transition, the polarization is 
lower than for a J = 1 — transition. We thus infer that 
the large polarizations (^ 10%) observed by us towards 
the Class I masers would require R ~ 100 F, intermediate 
magnetic fields (Ob/F ~ 100) and propagation preferen- 
tially perpendicularly to the magnetic field (cf. Fig. 3C;rf 
of NW90). A stronger polarization can only be achieved 
by considering bidirectional propagation. According to 
NW90, for J = 2 — 1 transitions, iZ ~ F, and propaga- 
tion perpendicularly to the magnetic field, ~ 10% linear 
polarization can be reached (for Ob > 10 F). 

In unsaturated masers, the largest fractional polariza- 
tions observed by us (pl ~ 30 — 40 %) can be reached by 
means of anisotropic pumping (Western & Watson, 1984, 
see however Elitzur, 1996, where maser polarization arises 
in the unsaturated growth phase). The radiation field that 
pumps CH3OH Class II masers is naturally anisotropic 
when it comes e.g. from a nearby star, or when the es- 
cape probability for the pump radiation is anisotropic. As 
remarked above, the latter situation may occur e.g. in a 
layer compressed by a shock. Anisotropic pumping pro- 
vides a substantial polarization: even for J^p > 1, and 
R ~ lO^'^F (unsaturated regime), and without magnetic 
fields, NW90 find, for axially symmetric pumping, a frac- 
tional linear polarization of 30 — 40%, with a position 
angle along and a maser propagation perpendicular to the 
symmetry axis of the anisotropic pumping. 

For Class I masers, thought to be coUisionally excited 
(Cragg et al., 1992), anisotropic pumping is possible if 
the collisions themselves are anisotropic. As suggested by 
NW90, such conditions may be encountered when the col- 
lision partners are charged particles undergoing ambipolar 
diffusion in a shock. However, it is uncertain whether the 
ionization fractions in maser regions are sufficient to ac- 
count for anisotropic coUisional pumping. 

In any case, we can also expect polarization due to 
an anisotropic loss rate F, for both maser classes. As 
mentioned above, the anisotropic absorption of IR pho- 
tons from the maser levels, populating torsionally excited 
states, will imply an unequal population of the magnetic 
substates of the maser levels. 
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4.4. Application to the observed mm CH^OH masers 

It is very likely that most of our CH3OH masers are un- 
saturated, given their weakness with respect to the lower 
frequency transitions of the same class. Unfortunately, de- 
termining the degree of maser saturation is notoriously 
difficult. In the following, we assume that the observed 
mm CH3OH masers are unsaturated. As discussed above, 
anisotropic pumping or anisotropic radiative losses are 
then a necessary condition to achieve the observed linear 
polarization. 

Although a quantitative analysis of CII3OH maser po- 
larization is still lacking, we can discuss our results in view 
of the models presented in the previous section. Figure 1 
shows the distribution of detections of fractional polar- 
ization among the sources listed in Tables 2 and 3. The 
highest linear polarizations (35-40 %) are reached in both 
classes. Although our sample is still too small to derive 
reliable statistics from it, we find that a large fractional 
linear polarization is regularly achieved in both classes of 
mm CH3OH masers. Whether anisotropic pumping and 
anisotropic radiative losses are needed or not to achieve 
such polarizations is a matter of current debate. 

4.5. Discussion of individual sources 

4.5.1. OMC-2 

The direction of the linear polarization of the 133 GHz 
Class I maser in OMC 2 is close to the outflow axis (Fig. 7) . 
The outflow direction is well defined both by CO obser- 
vations (Williams et al., 2003) and by the alignment of 
the 44 GHz maser spots (Kogan & Slysh, 1998). Since a 
shock at the "working surface" of the outflow will be too 
fast and destroy molecules, an oblique shock, propagating 
away from the outflow axis, explains both the required en- 
hanced gas-phase CH3OH abundance, and the asymmetry 
for an anisotropic loss rate (i.e. escape probability for IR 
photons). In this scenario, the magnetic field tends to line 
up roughly along the shock front, and only a maser prop- 
agation parallel to the shock plane will provide a strong 
maser gain. The magnetic field and the maser propagation 
can then be almost perpendicular, because this outflow is 
not viewed pole-on. Such a model provides an explana- 
tion for both the observed fractional linear polarization 
{Pl = 7% despite R ^ T) and its orientation parallel to 
Bsky and hence to the outflow axis (cf. NW90, their Fig. 
11, case (2)). 

4.5.2. DR21-W 

The case of DR21-W is less clear. Slysh et al. (1997) 
observed that the 84 GHz and 133 GHz lines have simi- 
lar spectral features. All the more it is astonishing that 
in DR21-W, both masers, similar in strength, line shape 
and fractional linear polarization have polarization angles 
that differ by 40°. The reason may simply be the vari- 
ability of unsaturated maser emission (the 84 GHz maser 
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Fig. 7. Schematic representation of methanol masers in OMC- 
2. Black dots: 44 GHz masers of Kogan & Slysh (1998). 
Contours: extent of the overlapping outflow system (in CO(l- 
0), from Williams et al., 2003). The dashed straight line is the 
outflow axis. Dotted contours: blue shifted lobes. Solid con- 
tour: red shifted lobe. The polarization angle of our 133 GHz 
measurement and its error are indicated by the shadowed line 
(whose length is the beamwidth, fwhm). The positions of the 
FIR sources are from Chini et al. (1997). 

was observed in February 2002, the one at 133 GHz in 
May 2002). If the polarization characteristics are constant 
in time, and if the spatial separation of the masers at 
84 GHz and 133 GHz is similar to that of the two 44 GHz 
components of Kogan & Slysh, this could be a hint at 
an inhomogeneous magnetic field and/or an anisotropic 
pump rate P or loss rate F with different symmetry axes 
for both masers (for a distance of 3 kpc, 2'.'8 correspond to 
8000 AU, cf. Plambeck & Menten, 1990, further references 
therein). However, the spectral profiles of the polarization 
angle of either source is flat, i.e. the magnetic fleld and/or 
P and F are homogeneous at least along each maser's axis. 



4.5.3. NGC7538 

In NGC 7538-Metl, the linear polarization is roughly per- 
pendicular to a suggestive bow-shock structure as traced 
by shock-excited H2 emission (Davis et al., 1998). This 
is schematically shown in Fig. 8. As for OMC-2, the po- 
larization angle proflle is rather flat, and a similar sce- 
nario may apply. The polarization angles of the 133 GHz 
Class I maser towards NGC 7538-Metl and of the 107 GHz 
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Class II maser towards NGC 7538-IRSl are very close to 
those of the A 870 /xm continuum emission from these posi- 
tions. Since the magnetic field direction in the dust emit- 
ting region and the maser region is likely to be the same, 
this may be taken as evidence for magnetic alignment of 
dust grains. However, in the unsaturated regime, it cannot 
be ruled out that the linear polarization is rather due to 
the amplification of polarized background radiation from 
dust. Then the fractional linear polarization of the maser 
radiation and its polarization angle are the same as for 
the continuum seed radiation (see e.g. Goldreich et al., 
1973, eq. 64, or Wallin & Watson, 1997, equations 17-19), 
whatever the explanation of the grain alignment may be. 
This is only valid if the magnetic substates are equally 
populated, i.e. the maser has to be weak, so that its own 
radiation does not yet affect the absorption coefficients. 
In order to fulfill this condition, we have to replace the in- 
equality i? ^ r by the more stringent i?/r ^ fle/AwD, 
where Awd is the Doppler linewidth. For a 10 mG field, 
and a mean maser linewidth of 0.5kms~^ (fwhm), we 
have rie/AwD ^ 10^^, since the CH3OH molecule is non- 
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Fig. 8. Schematic representation of the mm methanol masers 
in NGC 7538, centered at the position of IRS 1. Contours: ex- 
tent of the outflow system (in CO (1-0), from Kameya et al., 
1989). Dotted contour: blue shifted lobe. Solid contour: red 
shifted lobe. The polarization angle of our 107 GHz & 133 GHz 
CH3OH maser measurements (filled circle resp. filled triangle) 
as well as that of the nearest A 870 /im continuum polarization 
measurements (open circles, Momose et al., 2001) are indicated 
by the straight lines (whose lengths indicate the beamwidth, 
fwhm). The grey-shaded area indicates the H2 emission sug- 
gesting a bow shock (Davis et al., 1998). 



paramagnetic. In turn, R << 10^^ F is only valid for 
rather weak masers, otherwise the polarization due to the 
seed radiation will be destroyed by the (still unsaturated) 
maser itself. An observable quantity which allows to dis- 
tinguish between both regimes is the line-to-continuum 
flux ratio. Since this ratio was only measurable with lim- 
ited accuracy and spatial resolution towards two sources 
(see section 3), the data presented here do not yet allow 
for a meaningful conclusion. Interestingly, Pestalozzi et al. 
(2004) successfully model the 12.2 and 6.7 GHz Glass II 
methanol maser towards NGC 7538-IRSl by high amplifi- 
cation of background radio continuum emission. 

4.6. Spectral profiles of the polarization angle 

Most of the polarized Class I masers have flat po- 
larization angle profiles. Only two sources (S 231 and 
W51-Met2) have an almost constant slope (Fig. 9), of 
— 72?8/kms~^ and 97?6/kms~^ , respectively. As shown 
by NW90, a slope in the polarization angle profile is pro- 
duced by a gradient in the magnetic field direction, for a 
constant loss rate F and pump rate P. Consequently, the 
spectral profile of the linear polarization has a maximum 
at the velocity where the angle between the magnetic field 
and the maser propagation direction is largest (clearly vis- 
ible in S231 at 95 GHz, Fig. 3). 

Keplerian disks with an azimuthal, frozcn-in magnetic 
field component following the velocity vectors (as mod- 
elled by Barvainis, 1984, for the SiO maser in OrionIRc2) 
are an attractive model providing an unconstrained expla- 
nation for both a gradient of the magnetic field direction, 
and for an anisotropic pump or loss rate. For the case of 
the 6.7 GHz Glass II maser in S231, Minier et al. (2000) 
propose a masing disk as explanation of the kinematics of 
the maser spots. Pestalozzi et al. (2004) explain the main 
component of the 12.2 and 6.7 GHz CH3OH maser emis- 
sion towards NGC 7538-IRSl by a differentially rotating 
disk. On the other hand, De Buizer (2003) casts doubt on 
the disk interpretation, finding that the shock-excited H2 
emission towards young stellar objects associated with lin- 
early distributed methanol masers is mostly aligned both 
with the associated outflows and with the methanol maser 
spots (which should be aligned perpendicularly to the out- 
flow in the disk scenario). After all, a slope in the polar- 
ization angle profile can also simply be due to a changing 
direction of the symmetry axis for an anisotropic pump or 
loss rate, without any need for a particular magnetic field 
structure. Only more accurate models and higher spatial 
resolution observations will allow to distinguish between 
those scenarios. 

4. 7. Circular polarization 

The tentative detection of circular polarization in two 
Class I methanol masers (Fig. 10) is interesting and the 
mechanism producing the effect is not clear. We note that 
the Zeeman components of the CH3OH maser are unre- 
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Fig. 9. Polarization angle profiles of Class I CH3OH masers. 
The ordinates are all at the same scale, to better compare the 

end-to-end difference of polarization angles. The flat profiles 
(bottom row) are clearly distinguished from those with a slope 
(top row). 



solved: As estimated in section 4.5.3, Ob/Awd ~ 10~^. 
Due to the unsaturated maser emission, the circular po- 
larization observed in the 133 GHz masers of M8E and 
DR21-W can only be due to non-Zeeman effects, a con- 
clusion also supported by the absence of the classical S- 
shaped signature of Stokes V. As estimated above, we have 
fie ^ R and F, and the symmetry axis for the molecu- 
lar quantum substates (" quantization axis" ) in the masing 
medium is aligned with the magnetic field. Circular polar- 
ization in this regime was investigated by Wiebe & Watson 
(1998), and results from a misalignment of the local quan- 
tization axis and the direction of linear polarization. Such 
a misalignment can be achieved by changes of Bgky along 
the maser propagation path. Towards Class II masers, we 
did not detect circular polarization. Our most stringent 
upper limit is pc < 2.5% (for NGC 7538-IRSl, 107 GHz 
maser). Likewise, Koo et al. (1988) get a negative result 
from their survey of the 12 GHz line. 

4.8. The velocity components of Class II masers 

Another remarkable result is that the different velocity 
components of a given Class II maser, although different 
in flux density, display a similar fractional polarization 
and polarization angle. This suggests either masers that 
are not too far from each other (e.g. 500 A.U. in the case 
of the components A & B of Cep A East, see Figure 11), or 
that the velocity of a single maser is redistributed along its 
propagation direction, splitting the maser line into differ- 
ent spectral components (e.g. in Fig. 11 the components 




V|„ [km s '] 

Fig. 10. Detections of circular polarization towards M8E (up- 
per) and DR21-W (lower), both in the 133 GHz CH3OH maser. 

The circular polarization pc is indicated by open triangles, 
and the linear polarization pl by filled dots (both on the same 
scale). The dashed lines indicate the level of zero polarization. 
For comparison, the Stokes 1 spectra are shown as grey-shaded 
histograms. The axes of M8E and DR21-W are not on the 
same scale. 



labelled A in Cep A East, that arc spatially and spectrally 
unresolved by Mehringer et al., 1997). In both cases, we 
can expect a rather homogeneous magnetic field structure 
and/or distribution of the anisotropy of the pumping ra- 
diation field. As for velocity redistribution, Nedoluha & 
Watson (1988) show that such a breakup into multiple. 
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Fig. 11. Overlay of the linear polarization (filled circles, loft 
ordinate) and the polarization angle (unfilled triangles, right 
ordinate) with the spectrum (grey-shaded) of the 157 GHz 
CH3OH maser towards CepA. The projected distance of the 
components A and B is indicated (derived by Mehringer et al. 
(1997). 

narrow lines is even possible in presence of a smooth veloc- 
ity gradient of the order of a few kms~^ along the maser. 
Despite these changes in velocity, the magnetic field and 
- perhaps more importantly - the direction of the pump- 
ing anisotropic radiation field remain constant. A constant 
anisotropy parameter is expected if the maser length is 
small compared to the distance from the source of pump 
radiation. 

5. Conclusions 

1. We find linear polarization in 10 out of 14 Class I 
sources, and in 3 out of 7 Class II sources. The high- 
est polarization (39.5%) is found in the Class I maser 
L379, the most significant one in M8E (pL/crms > 25, 
both at 133 GHz). Both Class I and Class II mm CH3OH 
masers reach substantial polarizations well above 10%. 
The mm CH3OH masers, given their weakness relative to 
those at lower frequencies, are unsaturated. The polariza- 
tion observed possibly needs, especially for the larger an- 
gular momenta, anisotropic pumping and an anisotropic 
escape probability leading to an anisotropic loss rate. The 
maximal polarization of Class I masers is close to that 
observed for H2O masers. None of our three polarized 
Class II sources reaches a polarization as high as that of 
SiO masers. 

2. For Class I masers, the polarization angle profiles across 
the maser lines are mostly flat; two sources (W51-Met 2 & 
S 231) show evidence of a linear slope, suggesting that the 



direction of the magnetic field and of the symmetry axis 
for anisotropic loss or pumping changes along the prop- 
agation axis of the maser. The more common flat polar- 
ization angle proflles could be attributed to masers orig- 
inating from shocks that are seen edge-on, providing the 
largest possible gain lengths and a highly anisotropic es- 
cape probability. 

3. The polarization angles of the Class I and II masers in 
NGC 7538 show a remarkable agreement with the results 
from submillimeter continuum polarization, the latter be- 
ing generated by magnetically aligned grains. We suggest 
that the Class I maser is associated with a bow shock in 
the outflow. The magnetic field in this maser is likely en- 
hanced by compression; due to the shock, the loss rate is 
anisotropic, yielding a 133 GHz maser with a relatively 
strong polarization of 24%. The 107 GHz Class II maser 
towards IRSl is weakly polarized. For both masers, the 
polarization angles are close to those of the submm con- 
tinuum towards the respective positions of the masers. 
This provides evidence for magnetic alignment of dust 
grains, the magnetic fields in the dust emitting and mas- 
ing volumes having the same direction. Another explana- 
tion is the amplification of the polarized continuum emis- 
sion, which provides the seed photons for this unsaturated 
maser. In several of the polarized Class I and Class II mm 
CH3OH masers shown here, the line-to-continuum fiux ra- 
tio is suflaciently small to make such a scenario plausible. 

4. The polarization direction of the prominent Class I 

maser in OMC-2 is close to the outflow axis. This observa- 
tion can be interpreted as a maser originating in a lateral 
C-type shock compressing the gas in a plane roughly par- 
allel to the outflow axis, making the loss rate anisotropic. 
The linear polarization - relatively strong for an unsatu- 
rated maser - indicates a maser propagation perpendicu- 
larly to the magnetic fleld. The flat spectral profllc of the 
polarization angle lends additional support to this shock 
scenario. 

5. All Class II masers where we detect polarization have 

several (at least two) velocity components with the same 
polarization characteristics. This is indicative of velocity 
redistribution along the propagation direction of a single 
maser with a velocity gradient, and/or of nearby, but un- 
resolved individual masers with a constant magnetic field 
direction and a constant direction for anisotropic pump- 
ing or radiative losses. 

6. Significant circular polarization of —7.1% and 3.5% 
was tentatively detected towards two Class I masers at 
133 GHz (M8E respectively DR21-W). We attribute this 
result to a non-Zeeman origin, due to the weakness of 
the Zeeman splitting with respect to the linewidth of the 
CH3OH emission (the molecule is non-paramagnetic), and 
because the maser emission is unsaturated. 
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